Fluctuation and Dissipation in Liquid Crystal Electroconvection 



O 

o 

(N 
(N 



U 



o 
o 
o 

> 
X 



Walter I. Goldburg^, Yadin Y. Goldschmidt^, and Hamid Kellay^ 
^Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, PA 15260 
Centre de Physique Moleculaire Optique et Hertzienne 
U. of Bordeaux,!, 351 Cours de la Liberation, 33405 Talence, France 
(February 8, 2008) 

In this experiment a steady state current is maintained through a hquid crystal thin film. When the 
applied voltage is increased through a threshold, a phase transition is observed into a convective state 
characterized by the chaotic motion of rolls. Above the threshold, an increase in power consumption 
is observed that is manifested by an increase in the mean conductivity. A sharp increase in the ratio 
of the power fluctuations to the mean power dissipated is observed above the transition. This ratio 
is compared to the predictions of the fluctuation theorem of Gallavotti and Cohen using an effective 
temperature associated with the rolls' chaotic motion. 

PACS numbers: 05.40.-a, 61.30.-v, 05.70.Ln, 47.52.+j 



The fluctuation-dissipation theorem relates the spon- 
taneous fluctuations in a system at temperature T with 
its response to an external driving force. An example 
might be a simple harmonic oscillator surrounded by air. 
The air molecules exponentially damp its periodic motion 
and also cause the oscillating mass to fluctuate about its 
equilibrium position [Q. 

The experiment described here is concerned, not with 
a system in thermal equilibrium but rather, one which is 
in a steady state: an electric potential V applied across 
a weakly conducting liquid crystal (LC), generates con- 
vective motion within it when V exceeds a critical value 
Vc The LC is in good contact with its surroundings, 
so that it's temperature is fixed at the ambient value T, 
as it dissipates power P at a mean rate (P) . Our interest 
centers on the dependence of (P) and its rms fluctua- 
tions CTp as a function of V or, equivalently, the control 
parameter e = ((T^VK^) - 1) 

Recently, Gallavotti and Cohen (GC) Q have gener- 
alized the fluctuation-dissipation theorem (FDT) to in- 
clude systems that are driven far from equilibrium, like 
the one studied here. While their work motivated the 
present experiment, we are unable to verify its central 
prediction concerning the probability density function 
7r(P). They showed that, under appropriate conditions, 
the system, assumed to be driven into a chaotic state, ex- 
periences such large fluctuations in P that sometimes P 
can become negative, i.e. the driven, dissipative system 
momentarily can send energy back into the power sup- 
ply that generates the chaotic fluctuations. The theory 
makes a firm prediction about the ratio, 7r(P)/7r(— P). 

The GC theory is formulated in terms of the entropy 
production rate St = St, which is related to the average 
power gained by the system during a time period r by 
St — PT/kTss- Here the steady state temperature kTss 
is equal to the mean kinetic energy per particle. The 
authors show that if the chaotic motion in the system 
satisfies certain conditions, then 

7r(s^)/7r(-ST-) = exp(TS7.). 



The theorem, which GC call the fluctuation theorem 
(FT), has been generalized by Kurchan Q to systems un- 
dergoing Langevin dynamics and by Lebowitz and Spohn 
to general Markov processes The FT reduces to the 
FDT in the limit of vanishing driving force . 

For a macroscopic system like ours, it was not possi- 
ble to achieve such large fluctuations in P that render 
the power negative. Thus we were barred from directly 
measuring the ratio 7r(P)/7r(— P) predicted by the GC 
theory. What we did find were unexpectedly large values 
of the dimensionless standard deviation ap / (P) when V 
exceeded Vc- Our results become consistent with pre- 
dictions of the FT of GC, only if: (a) The temperature 
is taken to be the steady state temperature Tgs associ- 
ated with the quasi-particles chaotic motion and not the 
ambient temperature T associated with the microscopic 
thermal motion and (b) the correlation length of the ve- 
locity fluctuations is taken to be much larger than the 
typical size of the rolls. 

The liquid crystal used here is methoxy benzylidene- 
butyl aniline (MBBA) , a nematic that undergoes a tran- 
sition out of the amorphous state, when driven by an ac 
voltage, which will be written as V{t) = \/2Fsin(27r/<). 
Typically / is in 100 Hz range, and the measurements 
reported here were made at exactly this frequency. The 
electrodes, across which this voltage is developed, are 
thin, transparent layers of indium-tin oxide, enabling 
their visual observation. When V exceeds the critical 
value Vc= 6.0 Volts in MBBA, the convective motion of 
the LC commences. The resulting convective motion is 
analogous to Rayleigh-Benard convection, but with the 
applied voltage V replacing the temperature difference 
across the sample [||. 

Most studies of electroconvection in LC's center on the 
study of the patterns generated by the LC rather than the 
conductive behavior itself. Observations through a mi- 
croscope reveal that above Vc, the convective rolls form a 
stationary pattern. At a very slightly higher voltage, the 
dislocations develop and, with a further voltage increase, 
commence to move. At even higher values of V/Vc, the 
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orientational fluctuations take on a chaotic or turbulent 
appearance |^ and the orientational domains fluctuate 
rapidly, though their motion is slow enough for the eye 
to follow. Thus this driven system, that appears to be in 
a chaotic steady state out of equilibrium, is a candidate 
for checking the applicability of the FT. 

In the present experiment, the thicknesses d of the 
MBBA samples were 37.5 and 50 /zm. These samples 
were square with the lateral dimensions L ranging from 
mm to cm. It was possible to make samples with L less 
than 1 mm, but their resistance was so high (hundreds 
of MTl) that measurements of the fluctuations in P could 
not be reliably made. 

The experimental arrangement is extremely simple; 
the driving voltage V(t), obtained from a signal gener- 
ator, is applied across the sample, of resistance i?, which 
is in series with a much smaller resistance r, which in 
our case was usually 100 KSl. The time-varying voltage 
Vr is pre-amplified and then fed into a lock-in amplifier 
tuned to the driving frequency /. The lock- in output is 
fed into dc amplifier (PAR 113) before it enters the A/D 
input of a computer, which records the detected signal. 
The current through the sample, i = v^/r, and hence 
the power P — iV fluctuates slowly compared to 1//, as 
we will see. The full spectrum of the fluctuations can be 
captured by setting the lock-in time constant at r = 0.1 

Whereas this experiment is mainly concerned with 
fluctuations in P about its mean value, this latter quan- 
tity itself exhibits an interesting dependence on the con- 
trol voltage V . Because r << R, one might expect that 
(P) = V'^/R, so that a plot of ^/{P) vs V would yield 
a straight line of slope 1/R, and indeed this is so. But a 
more sensitive way of displaying the dependence of (P) 
on V is to plot the sample conductance G = (P) /V'^ as 
a function of V. 
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FIG. 1. Conductance G = (P) of MBBA sample in 

inverse Ohms. The control variable is the rms voltage, of 
frequency / = 100 Hz. The sample is a square of dimensions 
L = 4 mm and with d = 50 fim. 

Figure 1 is a graph of this type. For V < Vc — 6 V, 
G is independent of V ^ and the sample behaves like a 
resistor, with i? = 1/G = l/[6.6 x lO^^] = 150 UD.. This 
resistor-like behavior was verified by replacing the sam- 



ple with a metal-oxide film resistor having R = 100 Mil. 
One sees from this figure that the convective motion gen- 
erates increased energy dissipation in the sample when V 
exceeds = 6 V. 

This appears to be the evidence that a phase transition 
to a chaotic state enhances the entropy generation rate 
s-T (see above). The time interval t in this experiment 
is the integration time of the lock-in, which is r = 0.1 s. 
Note that the vertical scale in Fig. 1 spans a very narrow 
range in G. It is not surprising then, that the transition 
at Vc was not detected in a similar measurement of the 
conductivity as a function of V in MBBA B . 
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FIG. 2. Current through sample in nanoamps with its 
mean value subtracted out and a linear slow drift removed 
as well. The strongly fluctuating signal was recorded at \^ = 
1.67 Vc. The weak signal is at 1/ = 4 vohs = 0.667 Vc The 
lateral dimensions of this square sample are L=2.2 mm, with 
d= 37.5 /im. 

Turning now to the fluctuations in P about its mean 
value. Fig. 2 shows a time trace of i{t) at two applied 
voltages: V = 4.0 volts (weakly fluctuating signal) and V 
= 10.0 V (large, slow fluctuations). The sample thickness 
here is d = 37.5 fim, and its lateral dimensions are 2.3 
mm and 2.15 mm. Over the two- minute interval of the 
measurements, the mean current drifts slightly, an effect 
probably caused by a slight drift in room temperature. 
The sense of the drift varies in sign for different measure- 
ments. To compensate for this drift, a term linear in i is 
subtracted out. 

We have also measured the probability density func- 
tion of the power dissipation P at various values of V. 
There was no systematic departure from a Gaussian form 
for this function 7r(P); its kurtosis was quite close to the 
gaussian value of 3. The accuracy of these measurements 
is limited by the slow drift in the lock-in output men- 
tioned above. It suffices then, to report the width ap of 
this this function. 

Figure 3 shows the non-dimensional standard deviation 
(STD) of the power fluctuations 
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op/{p)^^{{p-{p)f)l{pf 

as a function of V . The sample is square with i = 4 mm 
and d = 50 /im. 



we find that the dimensionless STD is given by 
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FIG. 3. Dimensionless standard deviation of the power 
fluctutations, ap/ (P) lis V. For this sample L = 4 mm and 
d= 50 /.im. 

To interpret the power fluctuation data in Fig. 3, we 
assume that the noise that appears at V < Vc is un- 
correlated with the much slower fluctuations of interest 
above Vc and, being of small magnitude, we can neglect 
its contribution to ap/ (P). We will also assume that 
above Vc the bulk of the power is dissipated as a result 
of the convective motion. 

In the voltage interval where ap/ (P) is large, the fluc- 
tuation time is significantly longer than outside that in- 
terval. This observation was difficult to quantify, even 
though it is apparent from looking at a record if i(t). For 
example, the fluctuation rate appears larger than 10 Hz 
(the inverse of the lockin time constant) at V — i V, 
whereas it is smaller than 0.2HzatF = 8V. Atl4V 
the fluctuation rate has increased to roughly 1 Hz. 

The FC theorem of Gallavotti and Cohen concerns the 
rate of entropy production in thermostatted systems 
averaged over an observation time r. The power dissi- 
pation can then be written as P = StPss, where Tss is a 
temperature characterizing the kinetic energy of the fluc- 
tuations in the chaotic or turbulent steady state. Mak- 
ing certain assumptions about the time-reversal invariant 
chaotic motion that is generated by the driving force, Co- 
hen and Gallavotti obtained the FT theorem as stated in 
the introduction. If we further assume that the power 
fluctuations about the mean are gaussian. 



7r^(P) cx exp[-(P 
as they appear to be, we flnd that 
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Comparing this to GC result 
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(1) 



To interpret our measurements in the convective re- 
gion, V > Vc= 6.0 V, we will assume that the convective 
noise is uncorrelated with the noise that appears at, say, 
V = A volts, which is well below Vc, and we subtract it 
out in a definition of the excess noise EiV): 

E{V) = [ap{V) - ap{V = W)]/[{P) - {P{V = AV))] 

To account for the measured maximum value of 
ap/ (P) (~ lO^-^), we replace it with E{V) and define 
a "convective temperature" T^s by 



{l/2)Mv, 
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(2) 
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= exp(rP/fcTss), 



where M is an appropriate mass of quasi-particles par- 
ticipating in the chaotic motion, and Vrms is the charac- 
teristic velocity of the local convective motion of the LC. 
We emphasize that Tss differs from the molecular tem- 
perature, which is very close to the ambient temperature 
(~ 300 K), since the chaotic degrees of freedom are not 
in thermal equilibrium with the microscopic degrees of 
freedom. 

By observing the trajectories of small particles con- 
vected by the motion of the LC [^,0, one finds that, in 
order of magnitude, v^yis — voe, where vo = 16 fiia/s. 
Our own admittedly crude measurements of the motion 
of dust particles give Vrms — 400 ^m/s at y = 9.6 V for 
the 4 mm sample used in Fig. |^. 

The diameter d of a roll is roughly the plate spacing, 
and in our samples their length was typically five times 
as large. For the sample of Fig. 3 this gives M in Eq. 

~ 10~^ Kg, since the specific gravity of MBBA is 
close to unity. Taking Vrms ~ 100 /im/s one obtains 
2kBTss = 10^^ J. Using r = 0.1 s and the measured 
value of (P) at y = 10 V which is about 7 x 10~''W, 
the rhs of Eq. (jl|) becomes approximately 10~^, which is 
smaller by two orders of magnitude than the measured 
ap / (P) . To accommodate for this difference, we argue 
below that the motion of the rolls may be correlated over 
scales much larger than the size of an individual roll. 

The transition to the convective state is a continuous 
one as seen in Fig. 1; it resembles a continuous magnetic 
phase transition in a system in thermodynamic equilib- 
rium or a supercritical Hopf bifurcation |]ll| . The "or- 
dered phase" occurs for V > Vc- Recently it has been 
claimed |12 that a finite turbulent system is critical at all 
Reynolds numbers. It thus resembles a two-dimensional 
XY model which is characterized by a diverging corre- 
lation length throughout the ordered phase. This is be- 
cause all length scales between the size of the system, 
where the driving occurs, to the microscopic scale, where 
dissipation occurs, are connected by the cascading pro- 
cess. Thus, we would expect the system to exhibit a large 
correlation length (of the order of the system size) above 
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the onset of convection. This may provide a possible ex- 
planation for the large magnitude of the fluctuations. If 
the rolls exhibit a collective correlated motion, then the 
relevant mass could be as large as the total mass of the 
liquid crystal. This will increase our estimate of f2 suf- 
ficiently to bring the measurement and the calculation 
of E(V) into approximate agreement. The correlation 
length associated with the collective motion of the rolls 
should be distinguished from the much smaller length 
that is associated with the director fluctuations. This 
state of the LC in the present experiment is of course not 
turbulent, because the Reynolds number Re associated 
with the roll size and Vrms is much less than unity, but 
the motion is nevertheless chaotic. Hence the present 
experiment is very different from that discussed in Ref. 
Ip^ , where the Reynolds number was large. 

To establish whether the numerator in Eq. (|l|) is in- 
deed proportional to the mass M of entire sample, we 
measured ap / (P) and (P) in two samples having dif- 
ferent lateral dimensions, L ~ 2 mm and 11 mm. Both 
measurements were made at the same voltage V = 10.1 
V, so that Vrms presumably remains the same. Accord- 
ing to the above model, ap/ (P) cx Lj ^JJP). The ratio 

R ^ {apj {P))2mml{'7p/ {P))llmra - 1/2. This is tO bc 

compared with the measured ratio, R = 2. To resolve 
this inconsistency, we must assume that the correlation 
length of the power fluctuations saturates at a scale of the 
order of a mm. With this assumption, the measured and 
calculated values of R come into approximate agreement. 

In summary, we have measured the steady-state power 
dissipation (P) in a fluid system (thin layer of liquid crys- 
tal) in the vicinity of a phase transition from the quies- 
cent to the convective state, where the fluid motion of 
the MBBA liquid crystal presumably becomes chaotic. 
To maintain this chaotic motion at voltages exceeding 
the critical value, the system extracts energy from the 
power supply at an increasing rate, implying an increase 
in its mean conductivity G{V). 

The dimensionless width of the fluctuations, ap/ (P), 
is much larger than that associated with a resistor at am- 
bient temperature T but goes through a maximum near 
V ~ 2Vc. The magnitude of ap/ (P) near its maximum 
value can be approximately understood in terms of the 
fluctuation-dissipation theorem and the fluctuation the- 
orem of Gallavotti and Cohen, provided that one inserts 
into Eq. a temperature kTss that is identified with the 
kinetic energy of the quasi-particles excited by the flow. 
If the mass of these excitations is taken to be that of the 
convective rolls in the system, the dimensionless width 
of the power fluctuations is somewhat smaller than the 
measured value. An even better agreement is achieved 
if one assumes that the motion of the rolls is correlated 
with a correlation length much larger than the size of 
the rolls observed with a microscope. It thus appears 
that the mass of the excitations in a measurement of the 
power fluctuations is an appreciable fraction of the mass 
of the liquid crystal. 



We are grateful for informative contacts with G. 
Gallavotti, E. G. C. Gohen, J-B Yoo, J. R. Cressman, 
X. L. Wu and V. Horvath. Without the help of J. T. 
Gleeson, this experiment could not have been carried out. 
This work is supported by the National Science Founda- 
tion under grants INT-9603361 and DMR-9701668 and 
by the U.S. Department of Energy (DOE), Grant No. 
DE-FG02-98ER45686. 

Note added: After this work was completed we re- 
ceived a preprint from J. T. Gleeson et al. |jl^ showing 
a transition similar to that depicted in Fig 1. 
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